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Through appropriate projections of an exact renormalization group equation, we study 
^ ■ fixed points, critical exponents and nontrivial renormalization group flows in scalar fleld 



theories in 2 < d < 4. The standard upper critical dimensions dk = -^ttij k = 2, 3,4, . . . 
appear naturally encoded in our formalism, and for dimensions smaller but very close to 
Q-i! dk our results match the e-expansion. Within the coupling constant subspace of mass and 
r-| ! quartic couplings and for any d, we flnd a gradient flow with two fixed points determined 
by a positive-definite metric and a c-function which is monotonically decreasing along the 
! flow. 
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1. Introduction 

Wilson's exact renormalization group [1] provides a functional differential equation 
which dictates the way short-distance physics gets integrated into a long-distance effective 
action. This equation and its analogues (Wegner- Houghton [2] and Polchinski[3]) are cer- 
tainly powerful but too complex to be of practical use. Different approximations and pro- 
jections have been devised to bring these equations to more workable, yet non-perturbative, 
settings. 

In this Letter we elaborate on a projection of the Wegner-Houghton exact renormal- 
ization group equation due to Hasenfratz and Hasenfratz[4]. The basic idea is to focus 
on the evolution of the zero mode of the scalar field, with the highest momentum modes 
being integrated out into a self-interaction term of the constant mode. This setting retains 
non-linearities from the original equation and is amenable to both analitical and numerical 
studies. Extending previous work along the same lines [5], we analyze the fixed points of 
scalar theories for arbitrary dimensions. We then address the renormalization flows be- 
tween fixed-points and, in particular, the critical exponents. This somewhat conventional 
piece of work is then completed with a new study on the irreversibility of the flows. It 
is an open question whether renormalization group flows are gradient. A theorem due to 
Zamolodchikov[6] ensures such a property in two dimensions but only inconclusive work 
has been done in higher dimensions (see, for instance, [7]). Here, with some guesswork 
and brute force we are able to show that at first non-trivial order in our approach the 
flow is indeed gradient and thus irreversible. It is determined by a positive-definite metric 
in coupling constant space and a c-function which is monotonically decrasing along the 
flow, and it connects a unique Gaussian fixed point to a unique Wilson flxed point for any 
2 < d < 4. 



2. Projection of the Wegner-Houghton Equation 

The basic flow equation we will use is that of Wegner and Houghton[2] , while our nota- 
tion follows that of [4] . The basic and fairly intuitive procedure of Wegner and Houghton 
is to consider a generic scalar field theory action S*, with an UV momentum cutoff Aq. 
Starting from that, one then performs in the path integral an integration only over the 
outermost infinitesimal momentum shell, with momenta e~*Ao < Q < Aq, and t small. An 
effective action will result for the unintegrated fields, now with a slightly smaller cutoff, 
A = e~*Ao. Momenta in this effective action are then rescaled so their range again becomes 
< g < Aq, and with that the fields themselves and the measure in the action will scale 
with their appropriate scaling dimensions. When all this is put together, a differential 
equation results describing how the effective action changes as this second cutoff is lowered 
and more and more momentum degrees of freedom are integrated out. We do not give the 
derivation of this exact renormalization group equation, since it is presented in detail in 
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[2], but instead only state the final result: 

dS 1 r' f, d^S dS dS / d^S ' 



dt 2tjq \ d(p{q)d(p{-q) d(j){q) d(j){-q) \d(i){q)d(i){-q) ^ 



(2.1) 



where the prime in the first integral above indicates integration only over the infinitesimal 
shell of momenta e~*Ao < Q < Aq, and the prime in the derivative indicates that it does 
not act on the 5-functions in dS/d(f){q). 

By projecting onto the constant mode 0(0) of (piq), the above exact renormalization 
group equation is considerably simplified, and becomes the fiow equation only for the 
effective potential of the theory (this, of course, also projects out some terms which could 
contribute to the flow of the effective potential, but the fact that we are still able to find a 
very rich structure in the ensuing flow motivates this projection). We furthermore use the 
approximation used in Ref. [4] , constraining the effective action to have no other derivative 
pieces than the canonical kinetic term, that is, in coordinate space: 



S= I d'^x {-{d^cpf + V{cP)} . (2.2) 
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This leads to the following flow equation for the effective potential: 

V{x, t) = — ln(l + V"{x, t)) + d- V{x, t) + (1 - - - r])xV'{x, t) + const. , (2.3) 

where Ad/2 = [{A7T)^/'^T(d/2)]~^ , the dot is a scale derivative d/dt, x is the constant mode 
0(0), and we again refer the reader to the derivation in Ref. [4]. In the approximation we 
are using, Eq. (2.2), we actually leave no room for a wavefunction renormalization, and 
this turns out to imply that i] = above. For greater ease of calculations, we will actually 
study the equation for /(x, t) = V'{x, t), trivially found from the above: 

with 1] already set to 0. We remark here that the constant Ad can be absorbed by a 
rescaling of x, thus disappearing from the equation above, a fact we will make use of later. 
This is a reflection of universality in Eq. (2.4), whereby the shape of /* will depend on 
Ad but the critical exponents will not. This is the starting point of our calculations. From 
here one can proceed either by investigating numerical solutions [4] or by analytical means. 
We choose the latter, where we will use the following polynomial approximation for f{x,t): 

M 

f{x,t)^Yl ^2m-l{t)x^'^-\ (2.5) 
m=l 

with M an arbitrary integer (and, naturally, better approximations will have larger M), 
and where only odd powers are chosen because we want the potential V to be reflection- 
symmetric. This approximation has been widely used in the past (cf.[5],[8]). 
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3. Fixed Points 

With Eqs.(2.4) and (2.5) as our starting point, our first objective is to determine the 
allowed fixed point solutions and their properties. This is easily done by substituting a 
polynomial fixed-point solution 



M 



2m-l 



(3.1) 



with finite but arbitrarily large M into the fixed-point equation, i.e., Eq.(2.4) with / = 0. 
A Taylor expansion in x then leads to a set of M nonlinear algebraic equations, of the 
form 



J- ~r C]^ 



w,{c,,c,, c,) = (4 - d)c, + - = 



(3.2) 



wm-i{c*i,4,...C2m-i) =0 

wm{c*i,cI,...C2m-i) = 0, 

which can always be solved exactly and recursively up to Wm-i, giving c*^,c'^, . . . C2j^_i as 
a function of c|. That is substituted in wmi which then becomes a polynomial of order M 
in with the form: 

WM = k{d) cliaoid) + ai{d)c\ + . . . + aM-i{d)cl^~^) = , (3.3) 

with ai{d) being polynomials of order M — 1 in d. (Note that = ( ^ c*yi = 0) is 
always a solution for any d. This is the Gaussian fixed point.) As an example, for M — 6, 
we find: 

^^"^^ " 155925A5 
ao{d) = {d-4:){d- 3){3d - 8){2d - 5){5d - 12) 
ai(d) = -699456 + 899960ci - 436386d^ + 95973^^ - 8651d^ + 150d^ 
a2(d) = -8763072 + 9018200d - 3349144^^ + 527126d^ - 30234d^ + SOOd^ 
as{d) = -31764096 + 27691776d - 8408924^^ + 1030750^^ - 43166d^ + 300d^ 
a4{d) = -42872880 + 32815292c/ - 8479920^^ + 85U09d^ - 28049d^ + 150d^ 
a5{d) = -19261320 + 13251980d - 2990962^^ + 254333^^ - 6903d^ + 30d^ 
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For any given M, a complicated phase space of solutions {d,cl) can be found and 
plotted numerically. We have done this up to M = 7, and the important aspects of these 
solutions can be summarized as follows: 

i) The first important feature is that, rather unexpectedly, ao{d) always factorizes 

into 

MM 

n - = n - ^) = (d - 4)(d - 3)(d - -) • • • . (3.5) 

m=2 m=2 

This means that at the upper critical dimensions d — dk, k — 1,2, . . ., cl — is actually 
a double solution to wm = 0, which indicates a branching of fixed-point solutions below 
these critical dimensions. This is in perfect agreement with the multicritical fixed-point 
solutions known to exist below these dimensions. 

ii) To further corroborate the above, an e-expansion of Eqs.(3.2) and (3.3) about any 
critical dimension (with M = k) also leads to the known e-expansion solution given by 
Hermite polynomials, that is, for d = dk — e, 



r{x) = KkeH2k-iix/Xk) + 0{^), Xk = \j^-zr^ : (3.6) 

where is a constant depending on d^ (for instance, for k = 2, d2 = ^, and K2 = ^/A2/72). 
Note that a simple e-expansion of Eq.(2.4) will lead to a linear equation and thus cannot 
furnish this constant Kk- At higher orders in e we expect our results not to agree with 
the standard e-expansion since the present approximation does not allow for wavefunction 
renormalization. 

in) When M is increased by 1, the first M — 1 equations in (3.2) remain unchanged. 
The solution to the last one seems, as far as we have investigated, to lead to convergence 
of previous solutions as M gets larger (see also [5]). 

iv) The lower the dimension, the less trustworthy is this approximation or, conversely, 
the larger is the M needed. Altogether, we find that, for any d, some solutions represent 
true fixed points while others are spurious. For lower dimensions, the number of true 
nontrivial fixed points increases, but so does the number of spurious solutions. 



4. Flows and Critical Exponents 

Having found particular fixed point solutions in some approximation M, we can now 
study how the renormalization fiow approaches these solutions by determining the critical 
exponents. To find them, we study small i-dependent departures from some fixed-point 
profile f*{x): 

f{x,t) = r{x)+g{x,t), (4.1) 
where again a polynomial ansatz is chosen for g{x,t): 

M 

g{x,t)^Yl ^2m-i{t)x''^-' . (4.2) 

m=l 
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When (4.1) is substituted in Eq.(2.4) and only linear terms in g are kept, the resulting 
equation is 



// j_ / 

!•*/^^ (^ _i_ f*i\29 



(1 + /*/)^ (i + /*0 

With the ansatz (4.2) and matching powers of x in (4.3), we then find: 

M 

5i^^ nij{c*,d)6j, 



(4.3) 



(4.4) 



where flij is an M x M matrix which depends on the input values c* and d. For higher 
M, the entries are rather long and unwieldy; here we present as an example 0^ for the 
M = 3 case: 



/ 



V 



4Ad(9c 



-5c*(i+cn) 



(1+C*)3 

-18Ad cg(9cg^-10c*(l+et)) 



6Arf 
1+c* 

["i a) (l+c*)2 
18Ad(9c;^-5c;(l+ct)) 





20Arf 
1+c* 



(6-2d)-^y 



(4.5) 



The critical exponents will be given by the characteristic frequencies of the Eq.(4.4), 
i.e. the eigenvalues of Q. For the fixed points continuously connected to the only nontrivial 
fixed point above d = 3, we have calculated the critical exponents numerically up to M = 7 
for dimensions between 2 and 4 in steps of 0.1. For 2.9 < d < 4 our results are plotted in 
Fig. 1. It is worthwhile noting that as d ^ 4, the critical exponents merge with the tower 
of canonical dimensions (2,0,-2,-4,...), which are precisely the critical exponents of 
the trivial Gaussian theory at d = 4 (i.e., the canonical dimensions of the (0^, </)^, </)^, . . .) 
couplings in d = 4). This is an indication in our setting of the existence of a unique 
(Gaussian) fixed point at d = 4. We furthermore note that for d = 3 our two leading 
exponents 

1 

iy = — = 0.656 

^1 (4.6) 

a;2 = -0.705 

match fairly well results gotten by other methods (in [4], = 0.687 and UJ2 = —0.595; field 
theory calculations, high temperature expansions and Monte Carlo methods[10] all yield 
ly = 0.630 ± 0.003 and -1.0 < U2 < -0.5). 

It is also possible to perform an e-expansion on Eq.(4.3) around a critical dimension 
d = dk- We have done this for d — d2 — 4, calculating the critical exponents analitically in 
e to C(e^) and for all operators, and the procedure is identical for all higher multicritical 
points. To simplify our calculations we make use of the universality of the critical exponents 
in Ad and set it to one. Then, we substitute Eq.(3.6) for /* in Eq.(4.3) (and here it is 
important to have the correct coefficient Kk in Eq.(3.6), as well as the 0{e^) terms), and 
make the following ansatz for g{x,t): 



g{x, t) = exp{(a;f ^ + eu)''/' + e^<jjf')t} {go{x) + egi{x) + e^g2{x)) , 



,(1) I .2, ,(2) 



(4.7) 
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where uf^ = 2(2 — £), £ integer > 1, is any one of the critical exponents at d = 4. At zeroth 
order in e, one finds go{x) ~ H2i-i{x). At first order in e, using the previous solution for 
go{x), one finds an equation of the form: 

^^^^^"^^ ^ " ^""l + ^^^^ ~ ^0 ^'^""^ 

(1) ,1 (2£-l)(2£-3)A , , ^^'^^ 
+ (2 + ^ j -)) 9oix) + ■■■ , 

where Tig is the Hermite operator which annihilates H2e-i{x) and the dots are other 
Hermite polynomials. Since go{x) on the r.h.s. is precisely this Hermite polynomial, its 

coefficient must vanish for the equation to be consistent. This determines w^^^ ■ At next 

(2) 

order, is determined in essentially the same way, that is, by canceling on the r.h.s. a 
zero mode of the Hermite operator which appears on the l.h.s.. The final result is: 

= 2(2 - - I (1 + (2^-1X2^-3) ^ + 26^^^^^ _ ^^^2^ _ ^ ^^^3^ 

where £ = 1, 2, 3, .. labels the exponents for the different operators. At order e this agrees 
exactly with standard field theory calculations in the e-expansion[10]. At next order, as 
announced above, our result differs slightly from the standard one due to the absence 
of wave function renormalization and the truncation of the exact renormalization group 
equation itself. By a somewhat more cumbersome calculation it is also possible to find the 
0{e) correction to the critical exponents analytically in k for a generic multicritical point. 
The answer is: 

Critical exponents only characterize the flow very close to a particular fixed point. 
Another option we have is to study the flow globally by substituting Eq.(2.5) directly into 
Eq.(2.4). Matching powers of a; in a Taylor expansion leads to coupled nonlinear flow 
equations for Cj(t) in the form: 

Ci=w^{c), i = l to M, (4.11) 

where the Wi{c) are given in Eq.(3.2). Arguably, a polynomial ansatz does introduce a 
perturbative element into the essentially nonperturbative nature of renormalization flows 
between distant fixed points, and our approximation very likely misses some features of the 
true fiow. However, we believe that, again, the sensible and rich structure that emerges 
does justify the simplification. 

We have solved the nonlinear fiow (4.11) numerically with M — 3 in d — 3: 

3 C3 



Cl = 2ci + 



27r2 1 + Cl 



9 Co 5 C5 ^ ^ 

. _ 27 c| 45 C3C5 

C5 — 



27r2(l + ci)3 27r2(l + ci)2 
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The {ci{t), cs(t)) subspace of that flow is shown in Fig. 2. We note there the presence 
of a Gaussian (at (0,0)) and a Wilson fixed point, and a unique trajectory leading from 
the former to the latter. To determine that this flow is gradient and permits a c- function 
description is the object of the next section. 



5. c-Function 

We now study some features of the geometry of the space of local interactions. If the 
beta functions of a theory can be written as a gradient in the space of coupling constants, 

no = -9-'§^ (5.1) 

where g'^^ is a positive-definite metric, we know that the set of renormalization flows be- 
comes irreversible [9]. In such a case, there exists a function C of the couplings which is 
monotonically decreasing along the flows: 

making their irreversibility apparent, so that recurrent behaviors such as limit cycles are 
forbidden. In two dimensions it is possible to prove that the fixed points of the fiow are 
the critical points of C and that the linearized RG generator in a neighborhood of a fixed 
point is symmetric with real eigenvalues (the critical exponents). 

The renormalization group fiows found in the previous section are all well-behaved. 
Therefore it becomes natural to ask whether these fiows are gradient, i.e., whether there 
exists a globally defined Riemannian metric gij and a non-singular potential C satisfying 
Eq. (5.1). The general solution for an arbitrary number of couplings M would be extremely 
difficult. However, we find that it is possible to treat the case M = 2, namely, the subspace 
of mass and quartic couplings. The beta functions corresponding to the two couplings Ci, C3 
are given in Eq. (4.12) (where we restrict to C5 = 0). Because of the positivity of C3 (03 is 
the coefficient of (f)"^ in V and is required to be positive for stability of the path integral) 
it is appropriate to make the following coupling constant reparametrization: 

ci — > = ci 

(5.3) 

C3 ^ A = 6Ad C3 . 
In these new variables the beta functions take the form 



2m^ + 



dt 2 (1 + m2) 

d X _ (A-d) 3 ^ ' ^ 

A — — - 



dt 2 " 4(l + m2)2 
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and the fixed points become 



Gaussian : 



(m^,AG) = (0,0) 



Wilson : (m^^Xw) = ( 



4-d ^24 (4 - d) 
10 -d' 10 -d 



(5.5) 



)• 



Note that the Wilson fixed point merges with the Gaussian one at d = 4, similarly to the 
situation in Sec. 4. Now, by trial and error and considerable guesswork, the following 
solution to Eq. (5.1) can be found: 



C(m^, A) has the expected properties of a c-function: i) it has a maximum at the Gaussian 
fixed point, ii) it has a saddle at the Wilson fixed point, and in) there is only one fiow 
connecting both points ( we have not normalized the c-funtion to one for the Gaussian 
fixed point as often done in the literature). Naturally, this description corresponds to 
our particular parametrization in terms of m and A, which implicitly carries a choice of 
subtraction point. The variation of C between fixed points is reparametrization invariant 
and its positivity amounts to physical irreversibility of the fiow. A contour plot of C for 
d = 3 is given in Fig. 3, which depicts the space of theories in the basis given by m and A 
as a hilly landscape. 

For the sake of completeness, let us comment that the first mention of irreversibility 
of the renormalization group fiow was spelled out in the context of perturbation theory 
by Wallace and Zhia [9]. Later, Zamolodchikov[6] proved a theorem in two dimensions, 
the c-theorem, which relates the irreversibility of the flows to the basic assumption of 
unitarity in the Hilbert space of the theory. Several authors [7] have subsequently come to 
the conclusion that a similar theorem holds in any dimension in perturbation theory. More 
generally, any expansion where the space of theories is reduced to a manifold in a space of 
couplings will accomodate a c-theorem. Our setting in this Letter does not clearly fall into 
this category, due to the appearance of rational functions of the couplings in Eq. (5.4), 
and the explicit construction of the c-function, though to first non-trivial order, might be 
of relevance. 

A systematic approach to the irreversibility of the renormalization group fiow in the 
projected Wegner-Houghton equation should rely upon a computation of Zamolodchikov's 
metric {i.e. all two-point correlators between composite operators in the theory). This 
will require an exact renormalization group equation for the generating functional equipped 
with a source for composite scalar fields. 
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C{m', A) = ^(1 + m')' - |(1 + m'f + ]x\l + m'f 

^ O 4r 



3 A^ 
16 (4 - d) 



(5.6) 



and 




(5.7) 
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Figure Captions 

Fig. 1. Critical exponents for 2.9 < d < A corresponding to the relevant, marginal 
and the first two irrelevant operators in the M = 7 approximation. 

Fig. 2. d = 3 Renormalization group flows projected on mass and quartic coupling 
subspace in the M — 3 approximation, ci is plotted on the x-axis and C3 on the y-axis. 

Fig. 3. c- function contour of Eq. (5.6). The Gaussian point is at the top of the hill 
(0, 0), whereas the Wilson point lies on the saddle (—1/7, -\/24/7). 
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